Lagrangian statistics and flow topology in forced two-dimensional turbulence 



B. Kadoch * 
M2P2-UMR 6181 CNRS & CMI, 
Universite d^Aix- Marseille, Marseille, France 

D. del-Castillo-Negrete 
Oak Ridge National Laboratory, 
Oak Ridge, Tennessee, USA 

W. J. T. Bos 

LMFA-UMR 5509 CNRS, Ecole Centrale de Lyon 
Universite Claude Bernard Lyon 1-INSA de Lyon, Ecully, France. 

K. Schneider 
M2P2-UMR 6181 CNRS & CMI, 
Universite d'Aix- Marseille, Marseille, France. 
(Dated: September 16, 2010) 

A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence 
is presented. The topology is characterized using the Weiss criterion that provides a simplified 
tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic 
(deformation dominated) , and intermediate (turbulent background) . The flow corresponds to forced 
two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains with 
non-slip boundary conditions. In the double periodic domain, the probability density function (pdf) 
of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the 
flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the 
elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach 
and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of 
residence time in the topologically different regions are computed using the Lagrangian Weiss field, 
i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, 
the pdfs of the residence time have self-similar algebraic decaying tails. On the other hand, in the 
intermediate regions the pdf has exponential decaying tails. The conditional (with respect to the 
flow topology) pdfs of the Lagrangian velocity exhibit Gaussian behavior in the periodic and in the 
bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the 
Lagrangian acceleration in forced turbulence show a comparable level of intermittency in the periodic 
and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in 
all cases, by self-similar power law behavior with a decay exponent of order —2. 

PACS numbers: 47.27.E-, 47.27.T-, 47.27.N- 



I. INTRODUCTION 



The emergence of coherent structures is a well- 
established phenomenon in two-dimensional turbulence 
[1]. These structures are characterized by localized re- 
gions of intense vorticity and have a direct impact on 
the flow dynamics and the transport [2 . In particular, 
vortices tend to trap tracers for long times, while strong 
shear can lead to large relative displacements of tracers 
[SI H] . The combined effect of vortex trapping and free 
streaming along shear flows is important because it can 
lead to anomalous diffusion [5 , 6 . Therefore, it is of sig- 
nificant theoretical and practical interest to understand 
the relation between coherent structures and transport. 
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The present paper focuses on this problem from a La- 
grangian perspective. 

By means of direct numerical simulation, we study 
transport of passive tracer particles in a turbulent flow. 
Due to the nontrivial spatio-temporal dynamics of the 
flow field, the particle trajectories are stochastic and 
one has to resort to statistical tools to characterize 
the Lagrangian dynamics. This leads to the prob- 
lem of understanding the role of coherent structures on 
the Lagrangian statistics of tracers. Although impor- 
tant progress has been made, this problem is not well- 
understood. In particular, most studies have focused on 
the statistics of the Lagrangian displacements [3[ 13 [9] . 
Although this problem is of great interest for the under- 
standing of dispersion and transport of passive scalars, 
there are other Lagrangian quantities of interest. For ex- 
ample, the study of the statistics of the curvature and 
the acceleration of Lagrangian trajectories in fluids has 
been the object of recent experimental [8^ and numerical 
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studies flDHH]. 

The study of the acceleration and the curvature along 
particle trajectories in turbulent flows is directly related 
to the problem of intermittency and non-Gaussian statis- 
tics. This is because the Lagrangian acceleration and cur- 
vature time series typically contain "bursty" phenomena 
which give rise to slowly decaying non-Gaussian tails in 
the corresponding pdfs. In Ref. [8 a high speed parti- 
cle tracking system was used to measure the curvature, 
K, = |an|ii^^, where is the normal acceleration and 
ul IS the Lagrangian velocity, of trajectories in a tur- 
bulent closed water flow between counter-rotating disks, 
and it was shown that the probability density function 
(pdf) of the curvature values exhibits power law decay 
of the form ~ k~^I'^ . The same result was obtained in 
three-dimensional direct numerical simulations [10 . The 
statistics of the Lagrangian curvature and acceleration 
was also studied numerically in decaying two-dimensional 
turbulence in Ref. [11 where it was shown that the La- 
grangian acceleration exhibited exponential tails and the 
pdf of the curvature decayed as ~ k~^-^^ . 

The curvature of fluid particle trajectories is directly 
related to the Lagrangian acceleration. This quantity has 
received considerable attention during the last decade, in 
particular in the framework of homogeneous and isotropic 
turbulence [T2 [12 • The effect of boundaries on the La- 
grangian statistics was addressed in [T4^ where it was 
shown that, in the case of decaying turbulence, the La- 
grangian acceleration is more intermittent in conflned do- 
mains than in periodic domains. 

In this paper we focus on the influence of the flow- 
topology on the Lagrangian dynamics. We investigate 
the statistics of the Lagrangian velocity, acceleration 
and curvature, but in contrast to previous works, we 
study the problem of forced turbulence in unbounded and 
bounded domains. Most importantly, going beyond pre- 
vious studies on the relation between the flow topology 
and Lagrangian statistics, we consider the conditional 
probability of the Lagrangian curvature and acceleration 
on the flow topology which we characterize using the La- 
grangian Weiss criterion. 

The remainder of the paper is organized as follows. 
Section II deflnes the forced, two-dimensional turbulence 
model in periodic and bounded circular domains, and 
discusses the numerical methods used in the integration 
of the model and in the integration of the Lagrangian 
orbits. Section III discusses the characterization of the 
flow topology in terms of the Weiss fleld. The residence 
time statistics are presented in Section IV. The condi- 
tional statistics (on the flow topology) of the Lagrangian 
velocity, acceleration and curvature are discussed in Sec- 
tion V. Section VI presents the conclusions. 




FIG. 1: Snapshots of vorticity field at t=25 in forced two- 
dimensional turbulence in a double periodic domain (top 
panel) and in a bounded circular domain (bottom panel) . Red 
denotes cyclones, cj > 0; Blue denotes denotes anti-cyclones, 
cj < 0; and Yellow denotes regions of vanishing vorticity. 



II. TURBULENCE MODEL 

The turbulence model is based on the two-dimensional 
Navier- Stokes equations written in dimensionless form 

^ + V • Vc^ - v^^uj -F^ = -iv X (xV) , (1) 

V • V = 0, (2) 

where, V = (ix, v) is the flow velocity, = V x V is the 
vorticity, v is the kinematic viscosity, and an external 
force. We consider two different types of domains: a dou- 
ble periodic domain (unbounded) and a circular domain 
with no-slip boundary conditions (bounded). The term 
on the right hand side of Eq.Q is a volume penalization 
term which takes into account the no-slip boundary con- 
ditions using X as mask function [151 [161 [18] . This term is 
absent in the periodic domain calculations. In the case of 
bounded domains, the mask function vanishes inside the 
fluid domain (^ = 0), and is equal to one outside the fluid 
domain (x = 1), where the no-slip boundary conditions 
are imposed. The permeability 77 is chosen sufliciently 
small for given v in order to insure the convergence 
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of the volume penalization method [18]. Here we use 
zy = 5 X 10""^ and r] = 10~^. For the periodic geometry, 
the Reynolds number is Rcp = S\fE jv = 2.5x10^, where 
S = 27T corresponds to the domain size and E = 1/2 (V^) 
is the turbulent kinetic energy. For the circular geome- 
try the Reynolds number is Rcc = 2r\^/u = 2.4 x 10^ 
where r = 2.8 is the radius of the circle. 

To obtain a statistically stationary flow we consider 
= Fr — fSi/j^ where Fr denotes a random isotropic 
stirring at small {k = 8) wave- numbers. The Rayleigh 
friction term, — is added to avoid the accumulation 
of energy at large scales due to the inverse energy cascade. 
For the periodic domain we choose /3 = 1. However, for 
the confined domain case we choose (3 = since the wall 
plays the role of an energy sink and there is no need for 
extra large-scale dissipation. In the confined domain, a 
time correlation is introduced into the forcing term. In 
the calculations presented here, a discrete Markov chain 
is used as in Ref. [17 , = aF{tn) + {I - a^^^'^F^ 

where a = j^^^^^"" and Tc = 10~^. This time correlation 
is smaller than the viscous timescale, and only serves to 
improve the numerical stability. 

The numerical scheme is based on a classical pseudo- 
spectral method, fully dealiased, and the resolution is 
N = 512^ grid points. The time-integration is of semi- 
implicit type with At = 1.10~^ [18]. The particle trajec- 
tories were calculated by interpolating the Eulerian quan- 
tities and by using a second order Runge-Kutta scheme 
for time integration. The Lagrangian acceleration is the 
sum of the pressure gradient, viscous diffusion and exter- 
nal forcing 

aL = -Vp^uV^u^ f , (3) 

where / is constructed by applying Biot-Savart opera- 
tor to F^. The Lagrangian statistics were computed by 
ensemble averaging over 10^ trajectories. Each trajec- 
tory was integrated for 5 x 10^ time-steps, which corre- 
spond to 800 r, where r = 1/V2Z = 0.0622 is the eddy 
turnover time and Z = l/2{uj'^) the enstrophy. Figure [l] 
shows snapshots of the vorticity field at a fixed time in 
the periodic and the bounded domains. The character- 
istic emergence of coherent vortices in two-dimensional 
turbulence can be observed. Note the relatively strong 
vorticity values at the non-slip boundary in the circular 
domain. 



III. FLOW TOPOLOGY 

Since the flow is two-dimensional and incompressible, 
V = i X VV^, where t) is the streamfunction. The 

Lagrangian orbits are then obtained from the solution of 
the Hamiltonian system 

dx di/j dy di/j 

dt " ~d^' dt^ d^' ^ ^ 

where ip plays the role of the Hamiltonian, and (x, y) are 
the canonically conjugate variables. By flow topology we 



mean the topology of the streamlines, which corresponds 
to the topology of the Hamiltonian. The study of the 
location and dynamics of the fixed points of the Hamil- 
tonian, which are the stagnation points of the flow, is 
one of the most natural topological consideration with 
direct impact to transport. Elliptic fixed points give rise 
to trapping regions and hyperbolic fixed points induce 
rapid mixing and stretching [19 . However, the nontriv- 
ial spatio-temporal dependence of fluid turbulence makes 
the identification and tracking of the elliptic and hyper- 
bolic regions rather difficult. Among the experimental 
works addressing this problem are Refs. ^O] |2T] which 
proposed a method to characterize the elliptic and hy- 
perbolic points of flows exhibiting spatio-temporal chaos 
based on the measured Lagrangian curvature of tracer- 
particle trajectories. These, as well as other experimental 
studies (see for example the recent review in [E]), have 
pointed out the need for a deeper understanding of the 
impact of the flow topology on transport and the need for 
novel diagnostics to understand the Lagrangian statistics. 

The characterization of the flow topology used here 
is based on the Okubo- Weiss criterion [22l [23] which is 
an approximate method of partitioning a flow field into 
topologically distinct regions. From the physics point of 
view the flow is divided into vorticity dominated regions, 
which correspond to elliptic regions, strong deformation 
regions which correspond to hyperbolic regions and inter- 
mediate regions. The partition is based on the relative 
value of 

Q = s^-uj\ (5) 

where uj = dxV — dyU is the vorticity and s'^ = sf -\- is 
the deformation where 

du dv dv du 

dx dy ' dx dy ' 

Using the local value of Q, the flow domain can be par- 
titioned in three disjoint regions: 

i) strongly elliptic for which Q < —Qo 

ii) strongly hyperbolic for which Q > Qo 

iii) intermediate regions for which — Qo < Q < Qo 

where Qo is the standard deviation of the values of Q, 
i.e., Qo = \/ (Q^) where (•) is the ensemble average. Al- 
though this way of partitioning the flow has limitations, 
since it is based on the linearization of the Navier-Stokes 
equation [24^, it is a conceptually simple tool that has 
provided key insights in the understanding of the dy- 
namics of turbulence and Lagrangian transport in fluids 
and plasmas [3 , 25 . 

Figure [2] shows a snapshot of the spatial distribution 
of the three- level, normalized Weiss field, Q = Q/Qo, in 
double periodic and bounded domains, corresponding to 
the flow field in Fig. [T] The values of Qo are 170 and 
300 for the periodic and confined domain, respectively. 
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FIG. 2: Snapshots of three levels Weiss field, Q — Q/Qo, at 
t=25 in forced two-dimensional turbulence in a double peri- 
odic domain (top panel) and in a bounded circular domain 
(bottom panel) . The flow fields correspond to those shown in 
Fig. [l] Blue denotes strongly elliptic regions with Q < — 1; 
Red denotes strongly hyperbolic regions with Q > 1; and 
Yellow denotes intermediate regions with — 1 < Q < 1. 



Consistent with the fact that elliptic regions are domi- 
nated by rotation, the negative regions are concentrated 
inside cyclones, cj > 0, and anti-cyclones, a; < 0, i.e. the 
coherent vortices. The circulation cells surrounding the 
vortices correspond to strongly hyperbolic regions which 
are dominated by deformation. The background turbu- 
lent field is characterized by small intermediate positive 
and negative values of the Weiss field. Consistent with 
the no-slip boundary conditions, the boundary in the cir- 
cular domain exhibits high values of \Q\. 

IV. RESIDENCE TIME 

The space-time dependent function Q{x^ t) in Eq. ([5| 
is the Eulerian definition of the Weiss field. Alterna- 
tively, given a trajectory {x{t),y{t)) obtained from the 
solution of Eq. Q, the Lagrangian Weiss field is defined 
as the time dependent function Ql(^) = Q{x{t),y{t),t)^ 
i.e., the value of Q along the particle trajectory. Fig- 
ure [s] shows the probability density functions (pdfs) of 
the Eulerian Weiss field and the Lagrangian Weiss field 



for periodic and bounded domains. The Eulerian pdfs are 
obtained from a histogram of the spatial values of Q at 
a fixed time, and the Lagrangian pdfs are obtained form 
the histogram of values along the trajectories of a large 
ensemble of particles. The relatively good agreement of 
both pdfs indicates a consistent statistical sampling of 
the Lagrangian initial conditions. As Figs, [l] and [2] show, 
there are typically more coherent vortices in the periodic 
domain than in the bounded domain. This explains the 
asymmetry of the pdf in Fig. [3] towards negative Q values 
in the periodic case. An interesting difference is also ob- 
served in the decay of the tails. In the bounded case, the 
pdf exhibits slowly decaying algebraic-type tails indica- 
tive of very large intermittent values of the Weiss field 
coming from the boundary. In contrast, in the periodic 
case the pdf decays exponentially. 
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FIG. 3: Probability density function of Lagrangian (red solid 
line) and Eulerian (green dashed line) Weiss field in two- 
dimensional forced turbulence in a periodic (top panel) and 
bounded (bottom panel) domain. 



Figure [4] shows typical passive tracer trajectories in the 
periodic and the bounded domains colored by the instan- 
taneous value of the three- level Lagrangian Weiss field. 
The blue sections of the trajectories exhibit the expected 
spiraling motion resulting from the vortex trapping in el- 
liptic regions {Q < —1) which move around the domain. 
The yellow section denotes the incursions of the particle 
throughout the turbulent background, and red denotes 
the sections of the orbit in the circulating cells surround- 
ing the vortices. One of the main thrusts of this paper 
is the use of the Lagrangian Weiss field as a method to 
characterize the statistics of passive tracers. 
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FIG. 4: Typical passive tracer orbits in two-dimensional 
forced turbulence color-coded by the instantaneous value of 
the Lagrangian Weiss field in a periodic (top panel) and 
bounded (bottom panel) domain. Blue denotes strongly el- 
liptic regions with Q = — 1; Red denotes strongly hyperbolic 
regions with Q = 1; and Yellow denotes intermediate regions 
with -1 < Q < 1. 



In the double periodic and bounded domains, it is ob- 
served that particles tend to spend relatively short times 
in the strongly hyperbolic regions. This is to be ex- 
pected since these regions are unstable from a dynamical 
systems perspective. On the other hand, the relatively 
long stay of particles in strongly elliptic regions results 
from the trapping properties of vorticity dominated re- 
gions. To quantify these ideas we show in Fig. |5] the pdfs 
of the residence time, r, in the strongly elliptic, P^(t), 
strongly hyperbolic, P^(r), and intermediate, P*(r) re- 
gions. These pdfs determine the probability that a given 
Lagrangian tracer stays in a region with the same value of 
the three- level normalized Weiss field for a given time r. 
The Gaussian fluctuations characteristic of the turbulent 
background result in the exponential decay of P*(r), con- 
sistent with the results in Ref. [26 for three-dimensional 
turbulence. A possible explanation is that the particle in 
the intermediate zone can be related to a Poisson process 
which is characterized by a small time-correlation and an 
exponential pdf. In contrast, P^{t) and P^{r) exhibit a 
self-similar behavior, corresponding to algebraic tails in 
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FIG. 5: Probability density functions of residence time for 
bounded (green) and periodic (red) domains. Top: pdf, 
P^(t), in strongly elliptic, Q < —1, regions. Middle: pdf, 
P^{r), in strongly hyperbolic, Q > 1, regions. Bottom: pdf, 
P^{r), in intermediate, — 1 < Q < 1, regions. 



the pdf. This is likely caused by the weak chaos, sticki- 
ness, and strong correlations characteristic of the vortex 
cores and the circulating cells surrounding the vortices. 
The strong dynamical instability of the circulating cells is 
responsible for the considerably larger decay exponent of 
P^(r) which is consistent with the relatively shorter red 
sections {Q > 1) of the particle orbits shown in Fig. [4] 
Because vortex trapping is a local phenomenon insen- 
sitive to boundaries, P^{r) exhibits the same decay in 
the periodic and the bounded domain. However, since 
particles in the turbulent background wander through- 
out the entire domain P*(r) shows some dependence on 
the boundary. 

The proposed use of the Lagrangian Weiss field in the 
definition of the residence time could be of value in the 
study of non-diffusive transport in fluids and plasmas. 
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Early efforts in the study of passive scalar transport in 
flows exhibiting turbulence and chaotic advection were 
largely based on diffusion models. In these models, it 
is implicitly assumed that the particle trajectories can 
be described as an uncorrelated, Gaussian, Markovian 
stochastic process. However, several analytical, numeri- 
cal, and experimental studies have pointed out the limi- 
tations of these statistical assumptions, and recent works 
have focused on the development of more general mod- 
els, see for example Ref. [27 and references therein. A 
particularly productive approach is based on the use of 
continuous time random walk models, see for example 
Ref. [28 and references therein. These models generalize 
the standard Brownian random walk by allowing more 
general, non-Gaussian (Levy) pdfs describing the par- 
ticles' jumps and, most importantly, by incorporating 
waiting time distribution functions to account for par- 
ticle trapping. The knowledge of the jump distribution 
and the waiting time distributions opens the possibility of 
constructing macroscopic models of non-diffusive trans- 
port using fractional diffusion operators. The characteri- 
zation of waiting time distributions in particle transport 
studies in structured flows with simple spatio-temporal 
dependence is relatively simple, see e.g. Ref. [5 . How- 
ever, in turbulent flows the trapping structures exhibit 
random motion themselves and can appear and disappear 
in the course of the numerical calculation, e.g. Ref. 
When this is the case it is difficult to adopt an objective 
definition of a trapping event, and the construction of 
the waiting time pdf becomes nontrivial. One possible 
approach is to use the pdf of residence time in strongly 
elliptic regions, P^(r), to construct the waiting time dis- 
tribution. This pdf provides an objective quantitative 
measure of the time a particle stays on a vortex even 
when, as shown in Fig. |4j the vortex moves. 



V. CONDITIONAL LAGRANGIAN STATISTICS 

In this section we study the dependence of the statis- 
tics of the Lagrangian velocity, acceleration, and curva- 
ture on the flow topology, using the conditional pdfs of 
these quantities on the values of the Weiss field. For 
each field (velocity, acceleration, and curvature) we con- 
struct three conditional pdfs determining the probability 
of a specific value provided the particle is in a strongly 
elliptic, strongly hyperbolic, or an intermediate region. 
For example, P^{aL\Q < —Qo)daL gives the probability 
that the Lagrangian acceleration has a value in the range 
{aL^ciL -\- daL) while the particle is in a strongly elliptic 
region. The definition of the conditional probabilities for 
the other fields (velocity and curvature) and for the other 
regions of the flow (intermediate and strongly hyperbolic) 
is straightforward. 

Figure [6] shows the conditional pdfs of the Lagrangian 
velocity in the x-direction in periodic and bounded do- 
mains. In the periodic case, and have Gaussian- 
type dependence, i.e., in log- normal scale they exhibit 
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FIG. 6: Conditional pdf of Lagrangian velocity in the x- 
direction with respect to the three-level Weiss field in the 
double periodic domain (top panel) and in the bounded cir- 
cular domain (bottom panel). 



parabolic profiles. On the other hand, has an 
exponential- type dependence, as shown by the linear de- 
caying tails in the log-normal scale. This indicates that 
in the intermediate region comprising the turbulent back- 
ground, the velocity is more intermittent than inside the 
vortices and in the straining regions surrounding them. 
However, as the bottom panel in Fig. [6] shows, this dis- 
tinction is less clear in the bounded domain. The corre- 
sponding pdfs in the ^-direction (not shown) have very 
similar behavior. 

The conditional pdfs of the Lagrangian acceleration in 
the x-direction are shown in Fig. [7j It is observed that 
the level of intermittency (characterized by the heavy 
tails) is comparable in the double periodic and bounded 
domains. In particular, in all cases the tails exhibit expo- 
nential decay. This is in contrast to what is observed in 
freely decaying two-dimensional turbulence [14 in which 
the circular domain shows a significant higher level of 
intermittency. The conditional pdf for the intermediate 
values of the Weiss field, P*(aL| — Qo < Q < Qo), shows 
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FIG. 7: Conditional pdf of Lagrangian acceleration in the 
x-direction with respect to the three-level Weiss field in the 
double periodic domain (top panel) and in the bounded cir- 
cular domain (bottom panel). 
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FIG. 8: Conditional pdf of curvature with respect to the 
three-level Weiss field in the double periodic domain (top 
panel) and in the bounded circular domain (bottom panel). 



a sharp peak around zero indicating a stronger intermit- 
tency of the acceleration in the intermediate region. 

Figure [8] shows the pdfs of the conditional pdf of the 
Lagrangian curvature. In all cases self-similar power law 
behavior is observed with a decay exponent of order —2. 
As discussed in Ref. [8 , this can be explained as fol- 
lows. In isotropic turbulence, the Lagrangian velocity 
ul is close to a Gaussian distribution, which implies that 
uj^ is distributed accordingly to a chi-squared distribu- 
tion and thus is distributed according to an inverse 
chi-squared distribution whose tails decay as (l/ii^^)"^. 
Since = \an\uL~'^, this implies that, under the isotropic 
turbulence and Gaussian assumptions, the tails of the pdf 
of should decay as a^"^. Note that in three dimensions 
the same argument yields the exponent —2.5 [8 . The 
slight departures from this simple scaling observed in the 
numerical is a consequence of the nontrivial statistics and 
the dependence of the curvature on a^. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have presented a study of the rela- 
tionship between Lagrangian statistics and flow topology 
which we characterized using the Weiss criterion. The 
Weiss criterion provides a simplified tool to partition the 
flow field into topologically different regions. By flow 
topology we mean the topology of the streamlines. Since 
the flow is assumed two-dimensional and incompressible 
the streamlines corresponds to the level sets of the Hamil- 
tonian describing the Lagrangian orbits. Most of the pre- 
vious studies have considered the Weiss criterion in the 
Eulerian context. Here we focused on the Lagrangian 
Weiss field. That is, we considered the Weiss field along 
individual particle trajectories rather than as a function 
of the spatial coordinates in the flow. The flow field corre- 
sponded to forced two-dimensional Navier- Stokes turbu- 
lence in double periodic and circularly bounded domains. 
The study of the statistics focused on the residence time 
and in the conditional probability density functions of 
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the Lagrangian velocity, acceleration, and curvature. 

The flow topology was partitioned into strongly elliptic 
regions for which Q < — Qo? strongly hyperbolic regions 
for which Q > Qo, and intermed iate r egions for which 
—Qo < Q < Qo^ where Qo = \/ (Q^), is the standard 
deviation of Q values. In terms of the three levels Weiss 
field, these regions correspond to(3<— 1, (5>1, and 
— 1 < Q < 1. The strongly elliptic regions correspond to 
the vortices, the strongly hyperbolic regions correspond 
to the large strain regions surrounding the vortices, and 
the intermediate regions to the turbulence background. 

A very good agreement was observed between the pdf 
of the Eulerian Weiss field and the pdf of the Lagrangian 
Weiss field, confirming statistical convergence. In the 
double periodic domain, the pdf of the Weiss field exhib- 
ited a negative skewness consistent with the fact that in 
periodic domains the fiow is dominated by coherent vor- 
tex structures, i.e., elliptic regions. On the other hand, in 
the circular domain, the elliptic and hyperbolic regions 
seem to be statistically similar. The pdf of the Weiss 
field in the bounded domain is more intermittent than in 
the periodic domain due to the large vorticity and defor- 
mation fiuctuations at the boundary resulting from the 
non-slip condition. 

To characterize statistically the time that particles 
spend in the topologically different regions of the space, 
we computed the pdfs of the residence time, r, in the 
strongly elliptic regions, P^(r), in the strongly hyper- 
bolic regions, P^(r), and in the intermediate regions, 
P*(t). The pdfs P^{r) and P^{r) exhibit algebraic de- 
caying tails, and the pdf P^ij) exhibits exponential de- 
caying tails. P^{t) provides an objective quantitative 
measures of the time particles stays on a vortex even 
though, the vortex can exhibit a complicated motion. 
This motivates the use of this pdf to construct waiting 



time pdfs in continuous time random walk descriptions of 
anomalous diffusion in turbulent systems with coherent 
trapping structures. 

To study the dependence of the Lagrangian statistics 
on the fiow topology, we investigated the conditional pdfs 
of the Lagrangian velocity, acceleration, and curvature 
on the values of the Lagrangian Weiss field. The con- 
ditional pdfs of the Lagrangian velocity typically have 
Gaussian behavior in the periodic and in the bounded 
domain, with the exception of the conditional pdf on 
the intermediate values of the Weiss field in the periodic 
domain that shows exponential decay. The conditional 
pdfs of the Lagrangian acceleration have a comparable 
level of intermittency in the periodic and the bounded 
domains. This is in contrast to the freely decaying tur- 
bulence case, for which the circular domain exhibits a 
significant higher level of intermmitency. Finally, the 
conditional pdfs of the Lagrangian curvature displayed 
in all cases self-similar power law behavior with a de- 
cay exponent of order —2, reminiscent of the Gaussian 
character of the single point velocity statistics. 
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